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Abstract 

Myths and legends across the world contain many stories of deluges and floods. Some of 
these have been attributed to tsunami events. Doggerland in the southern North Sea is a sub¬ 
merged landscape and is thought to have been heavily affected by a tsunami such that it was 
abandoned by Mesolithic human populations at the time of the event. The tsunami was gen¬ 
erated by the Storegga submarine landslide off the Norwegian coast which failed around 8150 
years ago. At this time there were also rapid changes in sea level associated with the deglacia¬ 
tion of the Laurentide ice sheet and drainage of its large proglacial lakes, with the largest sea 
level jumps occurring just prior to the Storegga event. The tsunami affected a large area of 
the North Atlantic and left sedimentary deposits across the region, from Greenland, through 
the Faroes, the UK, Norway and Denmark. From these sediments, run-up heights of up to 20 
metres have been estimated in the Shetland Isles and several metres on mainland Scotland. 
However, sediments are not preserved everywhere and so reconstructing how the tsunami prop¬ 
agated across the North Atlantic before inundating the landscape must be performed using 
numerical models. These models can also be used to recreate the tsunami interactions with 
now submerged landscapes, such as Doggerland. Here, the Storegga submarine slide is simu¬ 
lated, generating a tsunami which is then propagated across the North Atlantic and used to 
reconstruct the inundation on the Shetlands, Moray Firth and Doggerland. The uncertainty 
in reconstructing palaeobathymetry and the Storegga slide itself results in lower inundation 
levels than the sediment deposits suggest. Despite these uncertainties, these results suggest 
Doggerland was not as severely affected as previous studies implied. It is suggested therefore 
that the abandonment of Doggerland was primarily caused by rapid sea level rise prior to the 
tsunami event. 


1 Introduction 

Global myths and ‘traditional oral tales’ contain many stories of great deluges or flooding events. A 
number of these have been attributed to tsunami inundation (Bryant et al. 2007, Heaton & Suavely 
1985, Antonopoulos 1992) or sea level rise (Nunn 2014, Giosan et al. 2009) and an increasing body 
of geological evidence is being discovered to test the validity of the legends (Bruins et al. 2008). In 
the southern North Sea, the now submerged island of Doggerland was still emergent at around 8 kyr 
and was occupied by Mesolithic communities (Ballin 2017). Due to rising sea levels Doggerland was 
eventually flooded completely (Coles 1998). However, studies have hypothesised that Doggerland 
was inundated by a large tsunami caused by the Storegga submarine landslide which occurred 
around 8150 years ago (Bondevik et al. 2012), with some estimates of 5m high waves impacting 
the island (Hill et al. 2014). If true, this would have been a catastrophic event, affecting the 
communities living there, and possibly causing the abandonment of the island (Weninger et al. 
2008). However, due to uncertainties in the relative timing of sea level rise events (Lawrence et al. 
2016) and the magnitude of the Storegga tsunami on Doggerland (Hill et al. 2014), it has been 
challenging to untangle the effect of these events on Doggerland. 
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The Storegga tsunami was caused by a large, 2400-3200 km^ submarine landslide (Haflidason 
et al. 2005) off the Norwegian coast (Fig. 1). The resulting wave spread across the Norwegian- 
Greenland sea, inundating the coastlines of Norway (Bondevik et al. 2003), the Faroe Islands 
(Grauert et al. 2001), Greenland (Wagner et al. 2007), the UK (Smith et al. 2004), and Denmark 
(Fruergaard et al. 2015). Wave heights of tens of metres are estimated to have inundated the 
Norwegian coastline and nearly 20m on the Shetland Isles (Bondevik et al. 2003). When a tsunami 
interacts with the coastline, the wave steepens and then travels inland. Depending on local geo- 
morphological characteristics, a wave can travel several kilometres inland. Typically, a tsunami is 
formed of several waves forming a wave-train. Each wave in this train will inundate inland depend¬ 
ing on the size of the wave. A number of variables can be used to measure the impact of a wave 
(Fig. 2). Tsunamis that have been observed are typically characterised by inundation distance or 
run-up height. Inundation distance can be derived using eye-witness accounts, debris lines or mud 
deposits from the wave. Run-up height is derived from watermarks left on buildings. Wave height 
is typically the height of the wave just offshore, before interaction of the tsunami with the seabed 
starts to become the dominant control on wave height. Both run-up and wave height are taken 
from mean local sea level. Once the wave has inundated it can take hours or days before the water 
recedes (Chini et al. 2012). This means with a typical period of minutes to hours a tsunami wave 
train can inundate land that is already flooded from the previous waves. 

The estimated wave heights or run-up heights of palaeotsunamis can be derived from sedimen¬ 
tary deposits left by the tsunami (e.g. Smith et al. (2004)). These deposits are often embedded 
within peat or other sedimentary layers and can be traced inland (Bondevik et al. 2003). Other 
deposits are found within lake sediments (Romundset & Bondevik 2011). From each deposit a 
minimum run-up can be inferred by calculating the palaeoaltitude of the site (Fig. 2) taking into 
account relative sea level changes from modern bathymetry/topography. The sediment deposit 
indicates the minimum inundation distance (and hence run-up height) as the grain size of the 
sediment decreases with increased distance from the shoreline as does the amount of sediment 
deposited. In the Japanese 2011 tsunami a 0.5 cm or greater sand layer was deposited for 57% 
to 100% of the inundation distance, depending on total inundation distance (Abe et al. 2012). 
The deposits formed by the Storegga tsunami in the UK are mostly found in the north-east of 
Scotland, including the Shetland Isles (Fig. 1). The Storegga deposits are generally a medium to 
fine sand layer, often pale in colour (e.g. Smith et al. 1999, Dawson & Smith 2000), but coarse 
sand to gravel in some locations (Long et al. 2016), found within peat layers (Smith et al. 2004). 
The deposits often contain peat rip-up clasts embedded within the sand deposit, indicating erosion 
(Dawson & Smith 2000). The sediment found in the tsunami deposits is clearly dependent on 
the local sediment characteristics, but also on the wave and local geomorphological characteristics 
(Takashimizu et al. 2012). 

The Storegga tsunami occurred at a time of rapid sea level rise (Clarke et al. 2004) - an 
estimated 1.4m over 500 years, which may have occurred as a series of pulses of decimetre to 
metre-scale (Lawrence et al. 2016). The primary cause of these sea level changes was the 8.2ka 
cold event whereby the rapid drainage of the Laurentide proglacial lakes caused disruption to the 
Atlantic Meridional Overturning circulation (Barber et al. 1999). These lakes were fed by the 
retreat of the Laurentide ice sheet, the largest of the northern hemisphere ice sheets, which in 
itself had a substantial effect of post-Glacial Maximum sea levels (Carlson et al. 2007). There 
is significant debate on exactly how and when the lakes drained; whether in a single event or as 
multiple events, and the precise timing of the events (Hijma & Cohen 2010, Lawrence et al. 2016). 
However, the 8.2ka event is generally thought to have been between 8740-8160 years ago (Lawrence 
et al. 2016), the most recent of which overlaps with the possible dates of the Storegga tsunami. 
These sea level variations had a major effect on Mesolithic communities including a period of 
migration (Turney & Brown 2007), which was undoubtedly related to rapid transgression in the 
region, with coastline moving landward over 30km at this time (Hijma & Cohen 2011). 

These rapid changes in sea level rise propagate uncertainties into the estimating of tsunami 
wave height from sedimentary deposits, as do uncertainties in the Glacial Isostatic Adjustment 
(GIA) models used, which can be of the order of metres at any given time slice (Kuchar et al. 
2012). In addition, when estimating the run-up height from sediment deposits several assumptions 
are made, including that there is no change in tidal range at the location (e.g. Smith et al. (2004)) 
which is generally not the case (Shennan & Horton 2002). Detailed knowledge of sedimentary 
changes that have occurred in the region over the last 8,000 years are lacking leading to further 
uncertainties, which include the preservation of coastal sediments (Forbes et al. 1995). 
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Figure 1: The Storegga final deposit (red) initiated off the coast of Norway and spread into the 
Norwegian-Greenland sea (main map). Sedimentary deposits in the UK are found in the Shetlands 
and northern Scotland and as far south as northern England (circles - main map and inset). The 
coastline at the time (black line) was different to today’s (filled landmass - UK only) due to changes 
in relative sea level. 


Although sedimentary deposits can be used to infer the tsunami wave height at select locations, 
they cannot be used to infer waves on now submerged landscapes. Numerical models can, however. 
Several modelling studies of the Storegga slide and tsunami have been carried out, which varied 
in resolution and numerical discretisation used, but all using very similar methods to generate the 
initial wave (Harbitz 1992, Bondevik et al. 2005, Hill et al. 2014). All models used a solid-block 
wave generator whereby a fixed shape is moved at a prescribed speed downslope. Deformation in 
the slide which can alter the initial waveform (Smith et al. 2016) is therefore not considered in such 
models. The two earliest models used modern bathymetry and differed in the shape of the wave 
generator (Harbitz 1992, Bondevik et al. 2005). The most recent modelling study includes the 
effect of relative sea level changes over the last 8,000 years to reconstruct the palaleobathymetry 
at the time improving the agreement between the modelled wave and the minimum wave heights 
derived from sedimentary deposits (Hill et al. 2014). The model of Hill et al. (2014) also used 
multi-scale resolution that varied spatially such that higher resolution was used along complex 
coastlines. However, several important effects, such as inundation, have not been considered in 
any of the models to date. All models had ‘cliff’coastlines where the boundary conditions used 
did not permit flow onto the shoreline. Thus flooding from waves approaching the shores is not 
modelled and the wave is reflected directly from the coastline back into the domain. For this 
reason all previous models compared the wave height some distance from this cliff-coastline with 
the run-up heights derived from sediments (Harbitz 1992, Bondevik et al. 2005, Hill et al. 2014). 
Assuming the slope and characteristics of the bathymetry are known, it is also possible to use an 
empirical relationship, such as that by Synolakis et al. (2009), to derive estimated run-up heights 
from these. However, Long et al. (2016) used such relationships for the most westerly Storegga 
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Figure 2: Terms used to describe the magnitude of a tsunami wave when inundating land. H is 
the tsunami wave height offshore, R is the run-up height, and D is the inundation distance. The 
relationship of these terms to the metrics inferred from tsunami sediment deposits is shown, as is 
the relationship between hypothetical current and palaeo-sea levels and tidal ranges. 


deposit found in the UK and found the estimated run-up height was lower (4.6 m) than the run-up 
height estimated from sediments (8-10 m depending on the relative sea level reconstruction used). 

Here, the multiscale modelling approach of Hill et al. (2014) is used to simulate the generation of 
the Storegga slide tsunami and propagate it towards the UK, including Doggerland. That simulated 
wave is then used to force the boundary of a second model, with increased spatial resolution, that 
is capable of modelling coastal inundation. Three separate inundation models are carried out on 
the Shetlands, the Moray Firth and Doggerland. The first two models are used to evaluate the 
agreement bewtten the coupled model to minimum run-up heights derived from sediment deposits. 
These coupled models constitute the first attempt to validate a full reconstruction of the Storegga 
tsunami from initiation to final inundation against inferred runup heights from sediment deposits. 
The reconstruction is then be used to recreate the wave that impacted Doggerland for which there 
are no sedimentary records. 

2 Methods 

2.1 Regional-scale model 

The regional-scale tsunami model used is that of Hill et al. (2014). Briefly, this model solves the 
three-dimensional non-hydrostatic Navier-Stokes equations on a spherical shell under a Bousineseq 
assumption in a rotating reference frame (Piggott et al. 2008). The equations are solved in a 
finite-element framework using a mixed discretisation of linear discontinuous Galerkin for velocity 
and a second-order continuous Galerkin formulation for pressure (Cotter et al. 2009). A combined 
pressure-free-surface boundary condition is applied on the upper surface (Funke et al. 2011) with 
coastlines and the seabed using a no-normal flow boundary condition with a velocity-dependent 
drag term. The simulation is a carried out on a spherical shell which is one element thick in the 
vertical. 

The tsunami is initiated using a dynamic solid-block model that has a prescribed shape and 
velocity-time profile. As the slide moves water is forced upwards at the leading edge and pulled 
downwards at the trailing edge, creating the wave. Shear along the top of the solid-block is also 
included. The tsunami is then propagated over a reconstructed palaeobathymetry derived from 
GEBCO modern bathymetry adjusted according to the GIA model of Bradley et al. (2011). For 
details of the set-up and parameters used, see Hill et al. (2014). 
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2.2 Local-scale inundation model 


To model the impacts of inundation, three local-scale models were constructed for the Shetland 
Isles, the Moray Firth and Doggerland, respectively (Fig. 3). The wave inundation model used here 
is a shallow water equation model using a finite element discretisation (Telemac2D). Telemac2D 
has been used to simulate inundation from tsunami waves previously (Horsburgh et al. 2008) and 
has a wetting-drying scheme that performs well at high resolution (Eric Jones & Davies 2006). 
It is capable of simulating the inundating flood wave celerity to a 3% accuracy in dam-break 
simulations (Hervouet 2000a), and has minimal numerical dispersion (Malcherek 2000). For each 
model a triangular, multi-scale mesh was constructed, based on modern bathymetric/topographic 
data which was then adjusted using the GIA model of Bradley et al. (2011) as per the regional-scale 
model. The mesh resolution is a simple function of water depth to roughly maintain a constant 
number of elements per wavelength and ensure the model complies with the CFL criterion. For each 
model, the boundaries of the computational domains and the mesh resolution distribution were 
constructed in QGIS (QGIS Development Team 2009) and qmesh (Avdis et al. 2017). qmesh was 
then used to convert both boundary geometry and mesh size metric to a suitable UTM projection 
and to construct a mesh, using GMSH (Geuzaine & Remade 2009). Similar to Horsburgh et al. 
(2008), a bed friction was represented as a Nikuradse roughness length of 0.01 m for low drag 
simulation, 0.05 for medium drag and 0.1 for high drag. These values are representative of clean 
gravel, light vegetation and heavy vegetation respectively (Hervouet 20006). A forcing boundary 
condition of each of these models was based on the free surface height and depth-averaged fluid 
velocity from the regional-scale model. Other boundaries were either solid or free-radiation (see 
below for details). For each model three different palaeobathymetric reconstructions were made, 
each based on the modern bathymetry but adjusted using Bradley et al. (2011) as in the regional 
model scenario plus an additional -5m, Om or +5m sea level change applied. These changes to the 
sea level reconstruction account for uncertainties in the GIA reconstructions. These were run for 
each of the three drag scenarios giving nine model runs in total for each of the three locations. 
Each model was carried out in a UTM projection. Using the available sedimentary data from 
the Shetlands and Moray Eirth regions, the model configuration that provides an inundation that 
most closely matches estimated wave heights can be derived and then used to assess if the Storegga 
tsunami was a catastrophic event on Doggerland. 




Eigure 3: The regional Storegga simulation mesh, with the three local inundation sub-domains 
highlighted. 


2.2.1 Shetland Isles model 

The inundation model for the Shetland Isles consisted of 30m resolution bathymetric data from 
SeaZone, complemented by 5m resolution topographic data from the Ordnance Survey, which was 
resampled to 30m resolution. Boundaries for this model were the 50m contour to reduce the 
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modelling of topography that was too high for the wave to inundate. Other boundaries are as 
shown in Figure 4. Mesh resolution varied from 60 m to 20 km, with all elements at bathymetry 
greater than 20 m having 60 m resolution. Bathymetry was used to alter mesh resolution according 
to the wave celerity {y/^) to ensure the wave could be adequately captured in shallow water. A 
total of 1,297,938 elements and 648,868 nodes constituted the mesh (Fig 4). 



Figure 4: Mesh and domain used in the local tsunami inundation simulations in the Shetlands. 
Inset shows close up and Yell, Unst and Fetlar. Thin white line is the palaeocoastline of the 
simulation with no additional sea level adjustments made. Note mesh extends inland depending 
on the topography in the region. Bathymetry data is ©Crown Copyright/SeaZone Solutions. 
All Rights Reserved. Licence No. 052006.001 31st July 2011. Not to be used for Navigation. 
Topographic data is ©Crown Copyright and Database Right (2017). Ordnance Survey (Digimap 
Licence). 


2.2.2 Moray Firth model 

As with the Shetland Isles model, the Moray Firth local-scale model was constructed from bathy¬ 
metric data from SeaZone and topographic data from the Ordnance Survey with an approximately 
30m resolution. Boundaries for this model consisted of the 40 m topographic contour and a smooth 
arc, from the north-western edge to the south-eastern edge, which acted as the forcing boundary 
for this domain. The minimum resolution used here was 45 m with a maximum resolution of 5 
km. Due to the fact that the majority of the domain was shallow water, no alteration of resolution 
based on bathymetry was carried out. Instead resolution varied smoothly from 45 m at the 0 m 
topographic contour to 5 km over a 10 km distance. All areas where the water depth was 3 m 
or less were set to the minimum resolution and the external boundary had increased resolution at 
1km which smoothly increased to 5km at 100 km distance. The mesh was constructed in a UTM 
30N projection space and had 503,696 nodes and 1,007,464 elements (Fig 5). 

2.2.3 Doggerland 

The area around Doggerland is on the boundary of the SeaZone data availability and the quality of 
data in the SeaZone dataset here is low. The bathymetry here was reconstructed from EMODnet 
data (Marine Information Service 2016). Nominal resolution of the data is around 188 m. The 
mesh is a simple box in a UTM 3IN projection space. The mesh had a minimum resolution of 1 
km and resolution varied linearly with depth to a maximum of 1.5 km. The final mesh consisted 
of 1,035,877 nodes and 2,071,752 elements (Fig 6). 
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Figure 5: Mesh and domain used in the local tsunami inundation simulations in the Moray Firth. 
Inset shows close up of the Inner Moray Firth and Cromarty Firth. Thin white line is the palaeo- 
coastline of the simulation with no additional sea level adjustments made. Note mesh extends in¬ 
land depending on the topography in the region. Bathymetry data is ©Crown Copyright/SeaZone 
Solutions. All Rights Reserved. Licence No. 052006.001 31st July 2011. Not to be used for Nav¬ 
igation. Topographic data is ©Crown Copyright and Database Right (2017). Ordnance Survey 
(Digimap Licence). 

2.3 Comparison to tsunami sediment deposits 

Given the uncertainties in the reconstruction of the palaeobathymetry, as well as in the drag 
coefficient and the Storegga tsunami itself, the simulations of the Shetland Isles and the Moray 
Firth were compared to the sedimentary record to assess the performance (Fig. 1). To carry out 
this assessment a simple metric, which was the percentage of sites where sediment deposits have 
been recorded for the Storegga tsunami were inundated in the simulation, was calculated across 
both the Shetlands and Moray Firth simulations. It is also possible from the model to derive 
maximum flooding depth (i.e. run-up height) across the simulated domain to visually assess how 
close the wave came to inundating each sediment deposit under the different assumptions made 
(i.e. drag coefficient and palaeobathymetric reconstruction). 

3 Results 

Results across the Moray Firth and Shetland Isles local-scale models show that the inundation 
height of the Storegga tsunami as calculated by the models is generally lower than that derived 
from sediment deposits, even when using a very low drag coefficient, as many sites, especially in 
the Moray Firth do not get inundated in the model reconstructions. Seven sites out of 21 where a 
sediment deposit has been found show a tsunami inundation in the low drag simulations. Increasing 
drag decreased the height of the tsunami, as would be expected, meaning fewer sites were inundated 
(a decrease to just two sites). However, note that given the number of Storegga sediment deposits 
which are found within peat horizons, which indicates the presence of vegetation, the “medium” 
drag models are perhaps the most realistic. 

Increasing sea level by 5m produced more sites that show tsunami inundation as would be 
expected, but still not all sites are inundated (15 of 21). It is also worth noting that the increase 
of sea level by 5m does not lead to a uniform increase in peak wave height showing the non-linear 
nature of tsunami inundation. Decreasing sea level by 5m results in fewer sites being inundated, 
as expected (three sites from 21), with no inundation in the Moray Firth. 

3.1 Shetland Isles 

The wave travels from the boundary, impacting the north-eastern coast of the Shetlands Isles 
(Unst) first, before travelling south, impacting Fetlar and Yell (see Video 1 in supplementary 
information). The first wave travels inland eastward, inundating palaeovalleys around the site of 
the Norwich deposit. There is also inundation from the north along a palaeovalley and widespread 
inundation on low-lying land on the eastern side of Unst. On the southern side of Unst, a small 
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Table 1: Maximum free surface height (MFS) and inundation coefficient (Ic = 0 for no inundation, 1 for inundated) for all nine modelled scenarios for each site 
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Figure 6: Mesh and domain used in the local tsunami inundation simulations of Doggerland. Thin 
white line is the palaeocoastline of the simulation with no additional sea level adjustments made. 
Bathymetry is derived from EMODnet data (Marine Information Service 2016) 


wave first floods from the north (the original wave refracted around the northern margin of Unst), 
followed by a much larger wave from the south, which is the original wave refracted around the 
southern coast. These two waves meet around the location of the two deposits at Snarravoe and 
Burragarth, which appears to not quite be inundated by this wave. This original wave then travels 
south and a period of slow regression of the wave occurs. Before the first wave has fully receded, 
the second wave arrives, inundating the eastern side of Unst, but not as far inland as the original 
wave. However, this wave is much larger on the northern edge and the timing is such that the wave 
refracts around the southern and northern margins of Unst almost simultaneously, inundating the 
land at sites around Snarravoe and Burragarth. Burragarth appears to be inundated primarily 
from the north and Snarravoe from the south, though the flow dynamics at this time are complex 
(Fig. 7). This wave is also much larger on Yell and inundates the palaeovalley that cuts across 
Yell east to west (see Video 1). 

On the Mainland of Shetland, it is the fourth arrival that generates the greatest inundation 
distance around Sullum Voe (sites 5 to 8). The first to third wave arrivals all cause some inundation, 
but it is the timing of waves approaching the embayment from the north and east that causes the 
largest flooding extent. At sites 9 and 10, four inundations occur with the first and fourth causing 
the maximum flooding extent (Video 2). 

Many of the sites experience a high level of flooding: up to 16m in the base sea level case. If the 
sea level is increased by 5m, then an additional three sites are inundated and the maximum flow 
depth increases to over 22m (Fig. 8). Note this is an increase of over 6m, despite a 5m increase 
in sea level. Nearby sites experienced an increase in flow depth of less than 5m between the two 
scenarios. Increasing the drag results in lower flow depths as expected and an increase in drag 
from low to medium is roughly equivalent to a 5m drop in sea level. A 5m sea level drop results 
in only three sites being inundated (those on Unst). 

The model here can also calculate bed shear stress. Here, a constant drag coefficient (0.0025) 
is applied to compute the bed shear stress in post-processing analysis. The results from site 1 
(Snarravoe) show peak bed shear stress capable of moving boulders (as in Japan 2011), but this 
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Figure 7: Bed shear stress at Snarravoe (site 1). Coloured bands indicate critical bed shear stress 
for different grain sizes. Solid line indicates magnitude of bed shear stress and arrows indicate 
direction (with north aligned with the y-axis and east with the x-axis - not scaled to magnitude). 


last for a very short time on the incoming wave (Fig. 7). The bed shear stress then rapidly drops 
on the outgoing wave and stays between 0.5 to 1 Nm“^, which is around the critical bed shear 
stress for medium to fine sand. At sites 2 and 3 the results are analogous, and again there is a 
sustained period of bed shear stress higher than the critical shear stress for fine sand. This is 
consistent with the deposits found there (Smith et al. 2004). 

3.2 Moray Firth 

The wave history within the Moray Firth is somewhat simpler than around the Shetland Isles. 
The first wave is small and causes little inundation. There is then significant draw-down from the 
negative leading edge of the much larger second wave. The inundation occurs over a longer time 
than the wave that impacted the Shetlands and subsequently the receding phase is also longer, with 
some areas still inundated at the end of the simulation, particularly near Dornoch Firth (Video 3, 
sites 13 and 14). A number of sites were not inundated in the Moray Firth simulations even with 
an additional sea level rise (Table 1). However, two of these were near the northern and eastern 
boundaries of the domain (11 and 21, respectively). Most of the hooding at the sites of interest is 
small, certainly less than Im fiow depth (Fig 9). 

As with the Shetlands, adding 5m to the reconstructed sea level means more sites are inundated 
(6 vs. 2), but again there is not an increase of 5m fiow depth at sites that inundate in both scenarios 
and in fact the increase is less than Im. In addition, many of the sites were close to hooding in the 
original simulation. A 5m drop in relative sea level shows no sites being hooded. Together these 
results show the non-linear nature of tsunami inundation predictions and the Moray Firth area 
requires either a larger wave or a slight increase in sea level in order to inundate the sedimentary 
deposit sites. 

3.3 Doggerland 

Based on the results from the sensitivity analysis from both the Moray Firth and Shetland simu¬ 
lations the low drag scenarios with either an additional 5m sea level increase or no additional sea 
level increase are the most plausible. The impact of the tsunami on Doggerland was therefore be 
assessed using these two scenarios. 

The wave impacts the northern coast of the islands, completely inundating the smaller chain of 
islands recreated to the north east of the main island (Fig. 10 and Video 4). The wave propagates 
up to 17 km inland in the north-east of the island and up to 21 km in the west of the island. 
Maximum fiood depth is around 7 m and a total of 2000 km^ was hooded (35 % of the land 
surface). If relative sea level was 5m higher, then maximum fiood depth increases to just over 9 m 
and the area hooded decreases to around 1200 km^ (60 % of the land surface). In both cases bed 
shear stress indicates that the tsunami would be capable of transporting cobble-sized grains (up to 
256 mm). None of the scenarios modelled show a catastrophic inundation of Doggerland where the 
entire island was hooded. Doggerland would have to be substantially smaller or the Storegga wave 
substantially larger for the entire island to be hooded. Based on the topographic reconstructions 
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Figure 8: Maximum flood depth for the low-drag, +5 sea level simulation, in the Shetlands showing 
sites 1-10. 


here, the maximum runup height to entirely flood the island would be 16.5 m in the standard sea 
level scenario, indicating an increase of 7.5 m in the simulated wave is required. 

4 Discussion 

The work here demonstrates that a simple-solid block slide model of the Storegga tsunami can 
recreate the fundamental aspects of the event. Most sites in the Shetland Isles and the Moray Firth 
are inundated by this model, but this does require both a low drag coefficient and a change to the 
estimated relative sea level at the time. By simulating the inundation of the Storegga tsunami, 
insights into how the wave inundated the land surface can be gained. In the Shetlands, Unst in 
particular, two of the three sites flooded occur where the is confluence of waves approaching from 
the north and south simultaneously. This may provide ideal conditions for depositing tsunami 
sediments due to a sudden drop in flow velocities (Fig. 7). The model shows this by rapid 
decreases in bed shear stress. Similar sites from the model could indicate locations where more 
sedimentary deposits could be found. The locations where tsunami sediments are found but were 
not inundated by the model are generally in areas of rapidly changing coastlines, e.g. barriers and 
spit development, or where very close to inundation in one or more scenarios. Increased detail in 
reconstrucitons of the palaeocostaline would undoubtedly improve the simulations presented here. 

The models also allow, for the first time, a recreation of the wave that inundated Doggerland. 
Whilst previous studies have focused on the catastrophic nature of the wave (Weninger et al. 
2008, Hill et al. 2014). Here it is shown that the wave likely flooded only a small portion of the 
island. The modern tidal range at Lowestoft (the nearest modern tidal station) is 1.94m. The 
wave inundated to a level similar to this. Depending on the timing of the wave the Storegga 
tsunami may not have flooded much more than a spring high tide. A spring high tide of 2m would 
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Figure 9: Maximum flood depth for the low drag, +5 sea level simulation, in the south-west corner 
of the Moray Firth showing sites 12-20. Sites 11 and 21 are off the map in the direction indicated 
by the arrows. 


inundate nearly 1500 km^ for the scenario considered with no additional sea level (25 % of the 
land surface) and around 900 km ^ for the scenario with an additional 5m of sea level rise (45 % 
of the land surface). In both cases the spring high tide is around 75% of the tsunami inundation. 
Other natural hazards could produce similar if not larger inundation, for example the 2013 storm 
surge in the region caused wave heights of up to 6m above datum (Spencer et al. 2014). Of course, 
this reconstruction is based on a simplistic reconstruction of Doggerland, which does not take into 
account changes in sedimentary regimes in the region in the last 8,000 years. Re-alignment of 
coastal regions during inundation show features such as barriers may be preserved, but in other 
cases these are transformed into thin veneers of beach sediments (Forbes et al. 1995). A full study 
of the sedimentological and geophysical data at Doggerbank is required to provide a more realistic 
recreation of Doggerland. 

The palaeobathymetric reconstruction is one of the main uncertainties in this study. Despite 
using a recent GIA model (Bradley et al. 2011), more recent models show different histories of the 
British-Irish Ice Sheet and hence show moderate differences in sea level histories around the UK 
(Kuchar et al. 2012). However, both GIA reconstructions don not replicate some recent sea level 
data from northern Scotland (Long et al. 2016). To compensate for this uncertainty, a sensitivity 
analysis of adding or subtracting an additional 5m of sea level change shows that the Shetland 
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Figure 10: Maximum flood depth on Doggerland using the zero sea level change (A) and +5m sea 
level scenario (B) and in both cases low drag. Green-brown colours show topography, blues show 
maximum flood depths. 


Isles reconstructions were less sensitive to this than the Moray Firth reconstructions, particularly 
if relative sea level was 5m lower. As with Doggerland, tidal ranges are not included in these 
simulations. Modern tidal range at Lerwick (Shetlands) is 1.73m. Again, the wave would have 
to arrive co-incident with spring high tide in order to compensate partially toward the increased 
relative sea level change required to better match sedimentary records. There are no modern tidal 
gauge stations in the Moray Firth, but Wick on the north-east corner of the Scottish mainland 
has a similar tidal range to the Shetland Isles. However, it is not certain that a wave arrives at 
the Shetlands Isles co-incident with spring high tide and can then arrive in the Moray Firth valley 
co-incident with the spring high tide there. A more detailed study including astronomical tides 
is required. These tidal changes are based on modern tidal ranges. Tides around the UK have 
changed in the past 8,000 years and are predicted to do so with current sea level rise (Felling & 
Green 2014). Estimates for tidal ranges at the time of Storegga are not available in the regions of 
interest, but could alter by up to 0.5 m (Shennan et al. 2000). A full reconstruction of tidal ranges 
at this time is required to reduce the uncertainty here. 

The method proposed here of comparing model output to sedimentary deposits by comparing 
if a site is inundated or not is perhaps the only viable method with current data. Comparing 
inundation between models and record with recent tsunami has been done using the observed 
inundation layer (Lpvholt et al. 2012). Similarly, for palaeotsunamis a comparison of the sand 
layer extent by either systematic coring or geophysical investigation would enable a much more 
thorough comparison between models and sedimentary data. It may also allow a recreation of 
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the wave offshore using numerical inversion techniques (Tang & Weiss 2015), which in turn could 
inform the initial tsunami source. 

An alternative interpretation of these results is that the wave generated by the solid-block model 
is too small as it impacts the UK. Rather than requiring that the wave impacts the UK at spring 
high tide, a larger wave would inundate further inland and flood more of the sites where tsunami 
deposits are found. However, a larger wave may also change the size of the wave as it impacted the 
Norwegian coast. Here, there are also maximum wave heights as a series of lakes have recorded the 
tsunami inundation. A higher lake did not record any tsunami deposit, so in places a maximum 
inundation height can be obtained (Romundset & Bondevik 2011). A deformable slide model may 
generate a larger wave (Smith 2016) and may increase the performance of this reconstruction. This 
remains computationally challenging due to the high resolution required around the slide and to 
model the coastlines in sufficient detail. 

Similar studies on the Thera (Santorini) eruption and tsunami in the Aegean sea show a complex 
picture of a large-scale event potentially being the cause of the collapse of the Minoan civilisation 
(Antonopoulos 1992). Modelling of the event does show the potential of a large tsunami inundating 
several coastlines (Novikova et al. 2011), though accurate dating of the sedimentary deposits may 
indicate the tsunami occurred a significant time before the collapse of the civilisation (Minoura 
et al. 2000). Likewise, in this study, timing is again key; with several rapid pulses of sea level rise in 
the ^ 500 years preceding the Storegga tsunami and the timing of Doggerland being abandoned. 
The modelling study here indicates that the ultimate cause of the flooding on Doggerland was 
likely sea level rise as the tsunami event was probably of comparable magnitude to storm events 
in the region. The pulses of sea level rise of decimetre to metre-scale in the 500 years prior to 
the Storegga event (Lawrence et al. 2016) likely inundated a substantial part of the island at the 
time and the tsunami wave was not large enough to be considered catastrophic. Studies similar to 
this and that of Novikova et al. (2011) could clarify the potential causes of other ancient flooding 
myths across the globe. 

5 Conclusions 

A recreation of the Storegga tsunami, including both relative sea level changes and inundation 
has been carried out around the UK for the first time. The results show good agreement with 
sedimentary deposits found for the Storegga tsunami in the Shetlands Isles and the Moray Firth. 
The model gives indications of bed shear stress changes, which can provide insights into the flow 
dynamics that can result in the deposition of sediment. The model can also be used to recreate the 
tsunami on submerged landscapes. Here, the model indicates the tsunami impacted Doggerland 
with runup heights of up to 9 m, causing up to 2000 km^ to be flooded (35 % of the land surface). 
However, this flooding may be comparable to tides in that area at the time, and almost certainly 
comparable to large storm surges in the region. These results therefore show that whilst a major 
event in Doggerland, it is probable it was not a catastrophic event and sea level rise was the most 
likely cause of the abandonment of Doggerland. 
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